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Quantum gates and simple quantum algorithms can be designed utilizing the diffraction phe- 
nomena of a photon within a multiplexed holographic element. The quantum eigenstates we use 
are the photon's linear momentum (LM) as measured by the number of waves of tilt across the 
aperture. Two properties of quantum computing within the circuit model make this approach at- 
tractive. First, any conditional measurement can be commuted in time with any unitary quantum 
gate - the timeless nature of quantum computing. Second, photon entanglement can be encoded as 
a superposition state of a single photon in a higher-dimensional state space afforded by LM. Our 
theoretical and numerical results indicate that OptiGrate's photo-thermal refractive (PTR) glass 
is an enabling technology. We will review our previous design of a quantum projection operator 
and give credence to this approach on a representative quantum gate grounded on coupled-mode 
theory and numerical simulations, all with parameters consistent with PTR glass. We discuss the 
strengths (high efficiencies, robustness to environment) and limitations (scalability, crosstalk) of 
this technology. While not scalable, the utility and robustness of such optical elements for broader 
quantum information processing applications can be substantial. 



I. PHOTONIC QUANTUM ALGORITHMS WITHIN A VOLUME HOLOGRAM 

It is the objective of this manuscript to argue that the construction of a small, lightweight, field deployable, 
inexpensive, low-dimensional, quantum computer is feasible using existing photonics technology. It is well recognized 
that photons offer great promise in quantum information processing (QIP) given their robustness to decoherence. 
However, it is this very resiliency that hinders their utility in quantum computing. Their weak coupling with atomic 
structures leads to substantial inefficiencies. While they may not find themselves as the central facet of a quantum 
CPU, they can, nevertheless, be an integrable part of the quantum information processing done by a quantum computer 
or in a quantum communication system where repetitive low-dimensional, fixed and robust quantum algorithmic tasks 
are required (e.g. quantum error correction, quantum relays, quantum memory/CPU buss, ...). 

Various approaches have been proposed for quantum computing. Perhaps the most familiar of these is the quantum 
circuit model (QCM) [1 . Here the challenge is to identify a suitable subset of a universal set of quantum gates that 
can faithfully represent a family of unitary operations on d qubits (d-partite), such that the the number of gates 
required scale polynomially in d, whilst the dimension grows exponentially as 2'^ . After the unitary evolution the 
output needs to be projected onto the computational basis and irreversibly measured. One realization of the QCM 
has been linear optical quantum computing (LOQC) and one-way quantum computing (OWQC), which are also 
known as cluster-state quantum computing [2j-jt5| . The essential feature of each of these approaches involves either a 
non- linear measuring process or a preparation of hyper-entangled input states (|/A^)), or both. Here a sequence of 
measurements are made to project the output state {\OUT)) onto one or another of the computational basis states, 
i.e onto a mutually unbiased basis (MUB). In this manuscript we will explore the utility of phototnic quantum circuit 
models. 

In this manuscript we suggest that in principle any quantum algorithm can be encoded within a single hologram, 
and that in practice, it has been found that one requires a Imm/N thickness of holographic material, where (N) is 
the dimension of the state space. This thickness will ensure high fidelity unitary transformations [Bj. Nevertheless, 
there are commercially available holographic media, e.g. OptiGrate's photo-thermal refractive (PTR) glass, that can 
support faithful holograms with a thickness approaching a few centimeters, thereby extending our considerations to 
10 to 20 dimensional quantum state spaces |7J. Every application of LOQC to quantum gates that we are aware of 
requires a cascade of interferometers which is highly unstable and impractical. The approach here will "lock" these 
interferometers within a tempered piece of glass that is resistant to environmental factors. 
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In this paper we will demonstrate that quantum gates and simple quantum algorithms can be efficiently simulated 
by the diffraction phenomena of a photon within a multiplexed holographic element. In Sec.|ll]we briefly outline the 
structure and function of a volume holographic grating that will be referred to throughout this manuscript. In Sec. [Till 
we describe the photonic quantum eigenstates we use. In particular we use the space spanned by the photon's linear 
momentum (LM) as measured by the number of waves of tilt across the aperture, i.e. the transverse component of the 
wave vector parallel to the face of the hologram. In Sec. |IV| we will briefly review our previous design of a quantum 
projection operator. In Sec.|V|we discuss the two properties of quantum computing within the QCM that make this 
approach attractive. First, any conditional measurement can be commuted in time with any unitary quantum gate - 
the "timeless nature of quantum computing." Second, photon entanglement can be encoded as a superposition state of 
a single photon in the higher- dimensional state spaces afforded by LM. Our theoretical and numerical results indicate 
that OptiGrate's photo-thermal refractive (PTR) glass is an enabling technology. Additionally, we give credence to 
the projection operator approach of Sec. IV by applying it to a quantum teleportation (QT) algorithm. We choose 
this algorithm because it exhibits both quantum entanglement and the timeless nature of quantum computing pQ [S] . 
We used coupled-mode theory and supportive numerical simulations in our research, all with parameters consistent 
with PTR glass "F, '5^. In Sec. VI we describe the construction of the familiar CNOT gate without the need for 
multiplexing. Our alternative approach will ordinarily require stacking more holographic gratings. In conclusion. 



(Sec. VII I we will discuss the strengths (high efficiencies, robustness to environment) and limitations (scalability, 
crosstalk) of this technology. The holographic approach presented here is not scalable, and with existing technology 
it is not reprogrammable; nevertheless, the utility and robustness of such optical elements for broader quantum 
information processing applications can be substantial. There should be many facets of a quantum computer other 
than the quantum CPU that could take advantage of this technology, including error correction, QKD relays, quantum 
memory management tasks and the like. 



II. VOLUME HOLOGRAPHIC GRATINGS 



Volume holography is used today for 2D image storage utilizing Z^A pixel s / jjLm? ^ which consumes only 1% of the 
theoretical volumetric storage density (l/A"^) [10 . This field was first introduced by Dennis Gabor in 1948 and it has 
been a popular research field ever since the development of the laser in 1960. Thicker holograms have more precise 
angular selectivity, i.e. its ability to differentiate the difference between two plane waves separated by a small angle, 
and under certain well-known conditions can achieve near perfect efficiency [llj . A hologram is considered a volume 
hologram if its thickness d ^ A^/A, where A is the characteristic period of the index of refraction of the grating, 
and A is the wavelength of the light. It is important for our purposes to emphasize that volume holography enables 
higher storage densities, and under suitable recording configurations can achieve near perfect efficiencies. This is not 
possible with thin holograms and planar gratings (spatial light modulators) where efficiencies can be at most 33 %. 

Briefly, the transmission volume holograms we consider here are formed when a "signal" wave, {r\S) — ^(r)e'*'-'') 
is directed into a holographic material and made to coherently interfere with an oblique "reference" plane wave, 
as illustrated in the left diagram of Fig. [T]for A^ = 3. In the figure the "signal" wave is a superposition of N plane 
waves, 

N N 
1=1 i=l 

where ai are pure phase angles. In this paper, we only consider planar reference waves, and the signal state as the 
superposition of plane waves. Ordinarily the signal waves will have variable phase and amplitude modulations. After 
the hologram is developed, and if we direct the identical signal wave, {r\S), into the hologram, then for a perfectly 
tuned hologram, the reference plane wave, {r\R), should emerge from the hologram as illustrated in the right diagram 
of Fig. [T] If it is not tuned, then other diffracted orders, e.g. modes parallel to the signal states, may emerge. 

How does one form a perfectly tuned hologram within the coupled mode theory? We have shown recently [Sj that 
near perfect efficiencies can be obtained if (1) the hologram thickness is tuned to its optimal thickness, (2) if the each 
of the signal's Fourier wave vectors have the same projection onto the normal to the hologram surface, i.e. they all lie 
on a cone with half angle 9s as shown in Fig. [2j and (3) each of the reference plane waves lie on their own distinct cone 
concentric with the first, with half angle Or and centered on the normal to the hologram face. We will also consider 
multiplex holograms wherein multiple independent exposures are made within the holographic material before it is 
developed. We demonstrated using coupled-mode theory that if the signal waves {Si\i_i^2....N form an orthogonal set 
under the L2 norm in the plane perpendicular to the waves propagation direction (z), i.e. 



{Si\Sj) = I S*{x,y)Sj{x,y)dxdy ^ 6ij, 



(2) 



FIG. 1: The left diagram shows a recording of a volume transmission grating by the coherent superposition of a plane reference 
wave \Ri} and a linear superposition of three signal waves 15") = e"*^ l^i) + e*"^ \i3S2) + e'"^ jSa) . On the right we show the 
function of the hologram. If the identically oriented signal wave 15") is sent into the hologram then the reference wave, \Ri}, 
will be reconstructed in the diffraction. The diffraction pattern will ordinarily consist of higher order diffracted modes parallel 
to the signal state. However, for a suitably tuned volume hologram perfect efficiency can be achieved, as shown in the right 
diagram [8] . This is why we constrained the signal wave components to a cone of half angle 6s centered on the normal to the 
hologram face. 

then perfect efficiency can be acliieved for each of the signals [5] . 

A volume multiplexed hologram that has achieved perfect efficiency (vifithin coupled- mode theory |12j ) under the 
"3+1" conditions outlined above provides a Hnear map between signal and reference modes. Physically it represents 
a projection (or redirection) operator or signal state sorter, |5] 

V = \Ri){Si\ + \R2){S2\ + ■■■ + \Rn){Sn\, (3) 

uniquely identifying each pair of signal and reference waves. 

Although the index of refraction within the hologram can be rather complicated, these devices are strictly linear 
optical components. Therefore, the diffraction patterns for a beam of photons will correspond exactly to the probability 
distribution for a single photon in the beam. For the remainder of this manuscript we will assume we are dealing 
with low number Fock states. We will describe in the next two sections (1) how we encode a d-partite quantum state 
onto a single photon that can be used in volume holography, and (2) how the unitary operation representing an entire 
quantum algorithm can be encoded in a multiplexed hologram. 

III. ENCODING A D-PARTITE STATE ON A SINGLE PHOTON 

The dimension, TV, of a quantum state space spanned by the direct product of d qubits grows exponentially, N = 2'^. 
Thus to represent a single d-partite state by a single photon requires vastly more quantum numbers than available 
to its two polarization degrees of freedom. However, the photon can be characterized by its extrinsic properties as 
characterized by its wavefront amplitude and phase. The potential of extending photon-based quantum information 
processing and quantum computing to higher dimensions was made possible in 1983 when Miller and Wheeler |13j 
described a fundamental quantum experiment utilizing a photon's orbital angular momentum (0AM), and when Allen 
et al [13] described Laguerre-Gaussian light beams that possessed a quantized orbital angular momentum (0AM) of 
Ih per photon. This opened up an arbitrarily high dimensional quantum space to a single photon [131 115) . Following 
these discoveries Mair et al. [531 US] unequivocally demonstrated the quantum nature of photon 0AM by showing 
that pairs of 0AM photons can be entangled using the non-linear optical process of parametric down conversion. 
Shortly thereafter, Molina- Terriza et al. [27j introduced a scheme to prepare photons in multidimensional vector 
states of 0AM commencing higher-dimensional QIP, with applications to quantum key distribution. Recently a 
practical method has been demonstrated to produce such MUB states using computer-generated holography with a 
single spatial hght modulator (SLM) [16]. 

While photons with specific values of 0AM have been emphasized recently in the literature, one can equally well 
utilize any other extrinsic set of orthogonal basis functions for higher-dimensional QIP, e.g. energy, linear momentum. 



FIG. 2: These Four LM states (arrows) that lie on a cone are used to generate a bipartite, or 4-dimensional signal state space. 
Here we chose the eigenstates (each a plane wave) in this 8-dimensional space with eigenvectors, \Si) = \9s,(j>i), i £ 1,2,3,4, 
here (f>i = 0, (f>2 ~ 7r/2, cps = tv and (f>4 = 37r/2. This can be generalized to higher dimensions, 2^*, by introducing d eigenstates 
around the cone, as long as the angular selectivity of the hologram is less than A(j> = 27r/2''. The wave vectors are defined in 
Eq.H 

angular momentum. While 0AM states respect azimuthal symmetry, transverse or linear momentum (LM) states 
respect rectilinear symmetry. Independent of the representation we use, the mutually unbiased (MUB) states will 
ordinarily be modulated in both amplitude and phase [16]. 

While the advantage of higher-dimensional QIP lies in its ability to increase bandwidth while simultaneously 
tolerating a higher bit error rate (BER) [T7], such transverse photon wave functions are more fragile to decoherence 
under propagation than the photon's spin wave function [I8j , and the divergence of the states in propagation may 
require larger apertures. However, since the spatial distances within a quantum computer are small this is not a 
concern here. This is not the case for quantum key distribution. Our preliminary results show that the efficiencies one 
can achieve with holographic elements exposed by LM plane waves and their superpositions (volume Bragg gratings) 
is substantially higher than that for photon wavefronts exhibiting 0AM [H] . Multiplexed volume Bragg gratings 
have been thoroughly addressed in the literature [7l [12j [191 EQ]. Since the efficiencies are potentially higher, and 
volume Bragg gratings are well understood, we therefore will restrict our attention in this paper to the orthonormal 
transverse photon states of LM. 

We concentrate in this section on the physics of transverse or LM states. We quantize the LM within the 2- 
dimensional plane of the face of the hologram, Fig. [2j We know that Bragg gratings can obtain perfect efficiencies as 
long as the component perpendicular to the face of the hologram {k±) for each of the photons' momentum are the 
same i.e. the plane wave states must lie on a cone. In addition the thickness of the hologram must be carefully chosen 

mm ■ 

If we are to represent a d-partite state by a single LM photon, we will require at least N = 2'^ LM states around 
the cone of half angle Og centered on the normal to the hologram's face. As a result, we can freely choose as our 
computational quantum state basis any 2'^ non-colinear plane waves. In this case, our integer quantum numbers 
will be the azimuthal angles, ^i, i S {1,2, ...2'*} locating the wave vector of the photon around the cone. Here 
we assume the hologram surface is in the x-y plane. For photons of wavenumber k these waves correspond to a 
transverse linear momentum = siTi{6a)hk = /ifcy, with components pf = sm{(j)i)hk]^]^, = cos (0i)/ifc|| and 
p| = cos (6)hk = hk± = constant; respectively. Here, fcy = k {\/ D) is the magnitude of the transverse component of 
the wave vector of a plane wave of wavelength A with one "wave of tilt" (X/D) across the aperture of breadth D. 

In the frame of the hologram, and in units where the speed of light is unity, the components of the 3-momentum, 
p = {p^ ,p^}, for each of our 2"* photons can be expressed in terms of their wavenumber, (fc), i.e., 

Pi^hki^h {kf,k\,kl}^hk {ain [0) con sin [9) sin {4>i), cos {9 s)} ■ (4) 

These TV = 2^^ plane wave states, {r\Si) — e*'^* '' define our computational basis for quantum information processing. 



MUB^{\Si),\S2),...,\Sn)} 



(5) 
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Each of these states represents a transverse Fourier mode of a photon; they are orthogonal under the norm 
{{Si\Sj) = 6ij) and span our 2''-dimensional state space. 

IV. REDIRECTING STATES IN A COMPUTATIONAL BASIS FROM REFERENCE TO STATE SPACE 

In Sec. [njwe outhned our recent progress in efhciently sorting a single photon with arbitrary complex wavefronts 
(Eq. |3])[8j 9J. And we showed that under certain conditions perfect efficiency) within coupled- mode theory) can be 
achieved. We showed that a multiplexed hologram could function as a d-partite version of the more familiar qubit 
polarizing beam splitter, i.e. a single optical element that can efficiently sort each of the d-partite states in the 
MUB, or computational basis (Eq. pi). As we argued, thick holographic gratings fortunately produce high diffraction 
efficiency in the first order [HI [121 [19]. If several predominant diffracted orders are required, as is the case for sorting, 
several independent fringe structures can coexist within the holographic material. Such multiplexed holograms have 
been used for multiple-beam splitters and recombiners |19| . and more recently for wide-angle beam steering [2 Ij. In 
a recent paper we described such a MUB-state sorter based on a multiplexed thick holographic element constructed 
from commercially available photo thermal refractive (PTR) glass [8] . Due to the unique properties of PTR glass the 
grating's thickness can approach several mm and be highly Bragg selective [55] . There is evidence that such sorters 
can be highly efficient, with Bragg efficiency > 95% [21]. Our recent simulations and empirical data on thick Bragg 
gratings indicate that they may provide an adequate solution to this critical sorting and beam steering problem for 
"single photon" QIP. 

While the signal state sorter, Eq.[3| can be constructed in a multiplexed volume hologram under the strict conditions 
outlined in Sec. [nj we could just as well construct a multiplexed volume Bragg transmission hologram, TZ, that 
recombines a beam. i.e. maps the reference waves onto their corresponding signal states. Mathematically that is the 
Hermitian conjug ate of Eq. (S) 

n^V^^ \Si){Ri\ + \S2){R2\ + ■■■ + \Sn){Rnl (6) 

This can be constructed physically using the same hologram that we used for V, but instead of directing the signal 
waves to the hologram to reconstruct the corresponding reference wave, we could direct a reference wave which will 
reconstruct the desired superposition of signal wave states, i.e. the signal beam, \S). For example, two identical 
multiplexed holograms (P) stacked back-to-back would represent the identity operator, X, on this state space. Any 
signal wave used to record the hologram emerges unchanged from the holographic sandwich; however, within the 
holographic material it will undergo a complicated evolution into the reference states, and back. Though, this 
identity gate may not be particularly useful in its own right, it does help illustrate our point. We will find in the next 
two sections that these holographic state redirection elements (Eq. [s]) can be useful in the construction of quantum 
computing gates. 

V. QUANTUM TELEPORTATION IN GLASS 

To best illustrate how a quantum algorithm can be encoded within a single hologram we will focus our attention on 
quantum teleportation (QT). This three qubit gate lives in an 8-dimensional state space. It exhibits all the essential 
features we need to address, i.e. entanglement and conditional measurements. The approach here will generalize, in 
principle, to any quantum circuit model. We will adhere as closely as we can to the notation found in Nielsen and 
Chuang [T]. 

The standard experimental realization of QT utilizes a pair of entangled qubits, |/3oo)j and an arbitrary qu5it state 
vector to be teleported, ip, where, 

IM = (100) + 111)) (7) 

V- = a|0)+/3|l) (8) 

The QT circuit consists of a CNOT gate, a Hadamard gate as well as two conditional measurements, one of which is 
fed into a Pauli X gate, the other into a Pauli Z gate as shown in Fig. |3| We exploit the " timelessness" of quantum 
computing by using the principle of deferred measurements and commute the two measurements in Fig. [3| with the 
quantum X and Z gates. We arrive at an equivalent circuit model for QT where all measurements are made in the 
end as illustrated in Fig. |4| We focus in this section on this equivalent tripartite QT circuit. Essentially QT amounts 
to (1) preforming a unitary transformation on the input tripartite state followed by, (2) two measurements. 



\OUT) ^Uqt\IN), 



(9) 
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FIG. 3: A quantum teleportation circuit |T]. Here the 8-dimensional tripartite input state |/A*') — ® |/3oo) is the direct 
product of the arbitrary qubit state, \ip) — l/\/2 (|0) + |1)) we wish to teleport, and a maximally entangled pair of qubit states 
|/3oo) ~ l/\/2(|00) + This circuit is composed of an initial CNOT operation, followed by a Hadamard gate, and two 

Pauli transformations {X and Z) conditioned on two measurements. 



where 

\IN) = IV) ® |/3oo). (10) 

This unitary transformation, Uqt, can be decomposed into a product of four sequential unitary transformations 
as illustrated by the dashed box shown in Fig. |4j Evolution proceeds from the left to the right. The input state 
first encounters (1) a CNOT gate on the first two qubits with unitary transformation, Ucnot, followed by (2) a 
Hadamard transformation on the first qubit with unitary operator, Uh, followed by (3) a Pauli X gate on the third 
qubit conditioned on the second qubit with unitary transformation, Ucx, and finally (4) a Pauli Z gate on the third 
qubit conditioned on the first qubit with unitary transformation, Ucz ■ We can therefore write the overall QT unitary 
transformation as a product of these four transformations, 

Uqt = Ucz Ucx Uh Ucnot (11) 



This 8-dimensional state space can be interpreted as a direct product of 3 qubit states (Eq. 10 1 with entanglement 
in the last two qubits (Eq. [7|. However, we can also represent the tripartite state as a pure state in an 8-dimensional 



state space without loss of generalization. In this way the entire input state (Eq. 10 ) can be represented by a single 
photon in an 8-dimensional LM eigenspace as described in Sec. |IV[ We just need to define a correspondence between 
a state in the computational basis and one of the LM eigenstates (Eq. [s]). We make this identification below. 

It is convenient for our analysis, and our application to holography, to represent the equivalent QT algorithm in 
matrix notation. Here we choose the standard representation of the qubits, 

\0)=(l) and |1)-('?V (12) 



^ J ' ' \^ 1 ^ 

The two input states can then be expressed as vectors in a 2 and 4 dimensional state space; respectively, 



1^) = ( ^ ) , (13) 



\M 



Therefore the input state can be expressed as a vector in an 8-dimensional state space, 




(14) 



\IN) = IV) (S> |/3oo) 
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FIG. 4: A equivalent quantum teleportation circuit pp. We use the "principle of deferred measurement" and commute the time 
ordering of the conditional measurements with the last two Pauli gates. This emphasizes the "timeless" feature of quantum 
information processing. The QT circuit thus becomes a single unitary transformation followed by two measurements and the 
transported state (top wire — >■ bottom wire). 



Following the definitions of Sec. |III[ we can freely choose the following computational basis vectors in this 8- 
dimensional state space and their correspondence to each of eight "signal" LM plane waves: 



iki -r 










Vo/ 



The quantum numbers of the "signal" state, 
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(16) 



\S,) = Ifc" 



(17) 



represents the projection of the linear momentum of the photon onto the face of the hologram, i.e. the transverse 
components of the LM. There would ordinarily be two components of this momentum, \Si) — |fc^,fc|), hence two 
quantum numbers. However, here each of these eight "signal" plane waves lie on a cone of half angle Og. Therefore, 
the first quantum number is constant for all states, 



— k'^ = constant, Vi. 



(18) 



We define the transverse quantum numbers by equally spacing the states around the cone's perimeter with 0i = 
and A(/) = 7r/4. This is a generalization to eight plane waves of the four wave configuration illustrated in Fig. [2] 
Therefore, one can say each state is described by two quantum numbers, 9 and 0; however, we hold 9 fixed so that our 
holograms can achieve perfect efficiency (within coupled mode theory) in modeling a given unitary transformation. 
This will be described below. 

The input state (Eq. 15) is then a linear superposition of four of the eight computational basis vectors; and 
correspondingly, can be represented by a linear superposition of four plane wave state vectors. 



1^4 



(19) 



This state has a wavefront with both amplitude and phase variations. In a similar vein we can define a corresponding 
set of eight "reference" plane waves that will be used to construct our holograms. Here we define them as we did for 
the "signal" waves, on a seperate concentric cone of half angle, 9r ^ 9s, 



(20) 
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Similarly, each of these reference waves live on a second cone of half angle 9r 7^ and are equally spaced, but offset 
with respect to the signal waves, around its perimeter with (f>i = n and A0 = 7r/4. It is important to note that the 
exact placement of the eigenstates around the cones and relative relationships between the reference waves and the 
signal waves are unimportant within the couple-mode theory. Any configuration can yield perfect efficiency. However, 
in practice it is preferable to make the angular separation between the signal and reference waves as large as possible, 
since this can lead to higher angular selectivity and higher efhciencies. This is not a concern in this paper. 

We can use these basis vectors to construct each of the four unitary transformations shown in th e da shed boxes of 
Fig.l4| thereby yielding a representation of the overall QT unitary transformation, Uqt, given by Eq.[ll] In particular 
we fmd, 
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(21) 



as well as, 
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(22) 



Therefore, in our computational basis the unitary transformation representing the QT can be obtained by substituting 
Eqs. [2lp2] into Eq.[TTJ 
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(23) 



This unitary transformation of the equivalent QT circuit shown in Fig. [4] accomplishes the desired transformation 

|1)-|0)\ ^ /|1) + |0) 



\OUT) ^Uqt\IN) 
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(24) 



thereby teleporting the state function from the top wire to the bottom wire. 

We describe how to encode the unitary transformation, Uqt into a single holographic element. By construction, 
the rows of Uqt form an orthonormal basis that spans the 8-dimensional state space. Each of these rows can be 



expressed as a superposition of the computational basis state vectors of Eq. 16 For this example, each row requires 
a superposition of exactly two of the eight computational basis state vectors. Therefore, each row can be represented 
as a linear superposition of the corresponding (Eq. |16[ ) plane wave states - two plane waves each for this particular 
quantum circuit. For example, if we consider the state vector representing each of the rows of the matrix in Eq. |23[ 
we can construct each of the the eight QT signal states, that we will use in the successive recordings in the 
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hologram, 
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(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



By, inspection we can reconfirm that these eight basis vectors form an orthonormal basis spanning the 8-dimensional 
state space. 

We construct a coordinate system where the z-axis is orthogonal to the face (the x-y plane) of the hologram. A 
faithful holographic representation of Uqt can be achieved by multiplexing eight recordings within the holographic 
material. The hologram is recorded so that each row of the unitary matrix, Uqt, is used to generate its own volume 
holographic grating. A qubit gate such as the Hadamard gate (H) is represented by a 2 x 2 matrix and would ordinarily 
require two recordings. A 2-qubit gate such as the controlled-not gate CNOT is represented by a 4 x 4 matrix and 
would ordinarily require four independent recordings (this will be discussed in the next section) . The qutrit QT circuit 
is represented by a 8 x 8 dimensional matrix (Eq. 23 ) and requires the following eight sets of independent recordings, 
between each of the "signal" states in Eqs. [25]|32| and its corresponding "reference" plane wave of Eq. [20]as follows: 

1. The first recording is made by the coherent superposition of "signal state (r|Si) and {r\Ri)] 

2. The second recording is made by the coherent superposition of "signal state {r\'S2) and (r|i?2); 



3. The eighth and final recording is made by the coherent superposition of "signal state (?^Sg) and {r\Rs) 

The recording geometry for the case of the fifth recording is illustrated in Fig. [5j 

Since (1) each of the eight signal states lie on a single cone with axis normal to the face of the hologram and form 
an orthonormal basis spanning the state space, and (2) each of the eight reference states lie on another concentric 
cone and span an isomorphic 8-dimensional state space, the oct-tuple multiplexed hologram (within coupled mode 
theory) faithfully represents the desired transformation, 

Uqt - ^ {\Rl){Sl\ + \Rl){S7\ + |i?2)(52| + \R2){Ss\ + |i?3)(54| + |i?3)(56| + |i?4)(53| + \R4){S5\- 
\R5){Sl\ + \R5){S7\ - \R6){S2\ + \R6){Ss\ - \R7){S4\ + \R7){Se\ - \R8){S3\ + \R8){S5\) ■ 

However, this operator maps quantum states in 'signal' space (inner cone of Fig. [5]) onto corresponding states in the 
reference space (outer cone of Fig. [5]) . In order to map this output back into the corresponding states in the signal we 



would need to ta ke t he diffracted output of the hologram described by Eq. 33 as input to a "redirection" as described 
by Eq. [g] in Sec. HI This can be constructed in a similar fashion by eight recordings, wherein the i*'* recording is 
made by the coherent superposition of the i*'' signal plane wave basis state {r\Si) — e''"'* '' with the i*'' reference plane 
wave basis state, {r\Ri) = e*'*'''', yielding, 

8 

^QT = ^|5z)(i?^|. (33) 

i=l 
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FIG. 5: Volume holographic design of the fifth of eight independent recordings used to generate the QT unitary transformation, 
Wqt, in PTR glass. 

Stacking the redirection hologram behind the hologram for the desired transformation will yield the desired QT 
operation. 

QT = 7^QT Uqt = ^ {\Si){S^\ + \S{){Sj\ + \S2){S2\ + \S2){S^\ + I ^3) (^4 1 +53) (^6 1 + |^4)(^3| + \Sa){S^\- 
\S^){Si\ + I^sX^tI - \S6){S2\ + \S&){Si\ - \S7){Si\ + \S7){S^\ - \Ss){S^\ + \Ss){S^\) 

VI. THE CNOT GATE: STACKING INSTEAD OF MULTIPLEXING 

One alternative to multiplexing is to make single recordings in each of many holograms and to stack the holograms. 
Here we provide a second illustrative example with a specific design of a quantum CNOT gate compatible with PTR 
glass. This gate is realized by stacking four holograms, and we will describe each below. The CNOT gate is a two 
qubit gate. Therefore the dimension of the state space is 4-dimensional. While this state space can be constructed as 
a product space of qubits by utilizing the polarization states of two correlated photons, it can also be represented by a 
single LM photon in a 4-dimensional state space. The CNOT gate can be constructed with a single photon. Following 
the arguments of the previous section, we freely choose four independent plane waves lying on the cone shown in Fig. [6] 
We associate to these independent transverse LM modes the four orthogonal quantum state vectors, \Sx), IS'2), 153) 
and 15*4). Any state vector, \ip), in this 4-dimensional state space can be written as a linear superposition of these 
states, 

1^) = a\Si) + I3\S2) + 71^3) + 5\Si), (34) 

with. 



(35) 
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FIG. 6: Volume holographic design of the 4-dimensional CNOT gate in PTR glass. The gate can be constructed by a stack 
of 4 LM gratings, or by a stack of two multiplexed gratings. 



Each of our basis states can be expressed in matrix notation, 



\Si) 



\So) = 



\S.) = 



, and 15*4) 



(36) 



In this computational basis the CNOT gate can be expressed by the following unitary transformation: 



CNOT = 



10 
10 
1 
10 



(37) 



If we let the z-axis be orthogonal to the face {x-y plane) of the hologram, the four volume holographic gratings are 
recorded by a suitable superposition of the set of four signal plane waves, 

{r\Si) = exp ^jfci • , (r|6) = exp {ik2 ■ r) , {r\S2) = exp {ikj, ■ , and {r\d) = exp (1/04 • r) , (38) 
and four reference waves, 

= exp [iKi ■ r) , (r|i?2) = exp {iK2 ■ r) , (j^^a) = exp («K3 • r) , and (r|i?4) = exp (1^4 • r) , (39) 
as shown in Fig. [6j 

The hologram is recorded so that each row of the unitary matrix of the CNOT gate is used to generate its own 
volume holographic grating. For a 2-qubit gate such as the CNOT gate we would ordinarily require four recordings; 
however, since the first two bits are just an identity matrix we need only two layers to transform the signal states into 
the desired reference states. In addition to one holographic recording per dimension of the state space, we also require 
the conjugate of the each grating (two in the case of the CNOT gate) in order to transform the diffracted reference 
waves from the reference waves back into the desired signal states. In particular, the CNOT-gate constructed from 
four holographic gratings stacked together as is shown in Fig. [6j 
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1. The first grating is recorded with the two coherent plane waves corresponding to states {Ss) and |i?4) . 

2. The second grating is recorded with the two coherent plane waves corresponding to states |54) and {Rs) ■ 

3. The third grating is recorded with the two coherent plane waves corresponding to states I-R4) and 15*4) . 

4. The fourth grating is recorded with the two coherent plane waves corresponding to states IR^) and jS'a) . 

The four gates will not diffract the first two signal states \Si) or |52). However, the first two gratings redirect the two 
signal states 15*3) and [54) into |i?4) and I-R3) respectively, in accordance with the Pauli X-gate, 

The first hologram is equivalent to the operator, 

Ui = \Si){Si\ + \S2){S2\ + \R4){S3\ + \Si){Sil (41) 

and the second hologram recorded with signal plane wave, (rlS^) and reference plane wave, (r^i?4), is equivalent to 
the operator, 

^2 = \Si){Si\ + \S2){S2\ + |i?4)(i?4| + \R3){S4\. (42) 

After the signal state , |/iV), passes through the first hologram the basis vectors change from {Si, S2, S3, S4} to the 
orthonormal basis, {Si, S2, R4, S4}. This explains why the third term in Eq. 42 is entirely within the reference space. 

While these two recordings could have been made in a single multiplexed hologram, we recover the same function 
by stacking the two together, thereby generating the CNOT operation, 

UcNOT^U2Ui - \Si){Si\ + \S2){S2\ + |i?3)(^4| + |i?4)(^3|. (43) 

However, the output of these two stacked holograms are the reference states, \Ri), 1^2), I-R3) and I-R4). In order 
to redirect these back to the proper signal states, we require the redirection operator similar to Eq. [6j This can be 
accomplished by recording a third hologram with the states |i?3) and jS's). The third hologram is equivalent to the 
operator, 

U3 = |^l)(^l| + \S2){S2\ + |5'3)(i?3| + \R4){Ra\- (44) 

Similarily, the fourth hologram is recorded with the states \R/i) and 15*4) and the fourth hologram is equivalent to 
operator, 

Ui = \Si){Si\ + \S2){S2\ + \S3){S3\ + |54)(i?4|. (45) 

Therefore, the combination of the four stacked volume holograms have the desired action - the CNOT gate. 

CNOT ^Ui-U3-U2 -Ui - \Si){Si\ + \S2){S2\ + |54)(53| + |53)(^4|. (46) 

One can apply these principles to design a universal set of quantum gates, as well as simple quantum algorithms such 
as QT. 

The advantage of stacking the holograms is that one can make the hologram thicker, thereby increasing the efficiency; 
however, achieving and maintaining the proper alignment should be more problematic. By multiplexing, we would 
need two holograms, each with two independent recordings in them. The first would be equivalent to the last two 
holograms in Fig.[6j while the second would be equivalent to the first two and would just redirect the reference beams 
into their corresponding signal states. The first two recordings are complementary to the second two - thus in some 
sense we are recording the "square root" of the CNOT gate. 



VII. DISCUSSION: ADVANTAGES AND LIMITATIONS 



Constructing simple quantum algorithms and quantum gates in volume holograms provides substantially greater 
optical stability than the equivalent optical bench realization. For LOQC, this stability is the overarching advantage. 
Often quantum operators, e.g. the simple projection operator given by Eq. [3j require a cascade of interferometers 
where the output of one interferometer is used as the input of the next [SI [55]. Therefore, as the dimension of 
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each state space increases, it becomes exceedingly hard to stabiHze and is simply impractical beyond two qubits. 
Other approaches, such as crossed thin gratings lack the efBciency needed for QIP. The device proposed here can 
potentially achieve this in a single piece of glass without the problem of misalignment. The technology presented 
here can potentially replace "fixed" optical components for a broad spectrum of classical and quantum photonics 
experiments. 

The primary limitation of volume holographic QIP is that it is not scalable. Experience shows that multiplexing 
requires approximately 1mm per recording of the state space to achieve high fidelity, and in QIP applications this 
scales exponentially with the number of qubits. But then again, we are not aware of a realistic scalable quantum 
computer to date. Secondly, the holograms discussed here are write-once holograms and cannot be erased. Therefore 
the algorithm is "fixed" into the holographic element. While there are re-recordable holographic media, none that 
we know of has the specifications to outperform PTR glass for the applications discussed in this manuscript. This is 
hardly a prescription for a quantum CPU; however, as mentioned in the paper, this technology might be integral to 
complete QIP systems where smaller d-partite operations are needed on a routine basis, e.g. a quantum memory bus, 
quantum error correction circuit, QKD relay system, etc... 

While we have extensively analyzed volume holographic gates and algorithms using coupled-mode theory, paraxial 
wave equation simulations and finite-difference time domain simulations, we have not analyzed the engineering partic- 
ulars of this device. For example: (1) how many independent writes of orthogonal states into a holographic element 
can be made in the PTR glass before cross-talk or saturation between the modes becomes a limiting factor? (2) Is 
it difficult to stack the holograms due to the enhanced angular selectivity of the volume holograms? And (3) what 
is the maximum number of recordings in a multiplexed PTR hologram that can be reasonably be achieved? In this 
sense, we are well along in understanding these devices from a theoretical prospective; however, we are at the very 
beginning experimentally. 
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